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1 Introduction 


For finite collective velocities where the adiabatic approximation for the single- 
particle motion is no longer valid, a qualitatively new feature of collective nuclear 
motion has been predicted : elastoplasticity [1,2,3,4]. This dynamical behaviour re- 
sults from a coherent coupling between collective and intrinsic degrees of freedom and 
subsequent equilibration by residual two-body collisions. Whithin a non-Markovian 
transport-theoretical approach [3] , the elastic response is described by scaling of di- 
abatic single-particle wave functions [5] according to the collective deformation, while 
two-body dissipation is accounted for by a relaxation ansatz. This dissipative diabatic 
dynamics (DDD) ascribes elastoplasticity to nuclei and establishes a link between time- 
dependent Hartree-Fock and Markovian transport theories of nucleus-nucleus collisions 
[6]. Isoscalar giant quadrupole vibrations of nuclei and mass diffusion in nucleus-nucleus 
collisions can be considered as well-established examples for the elastic and plastic lim- 
its of elastoplasticity. The first numerical results obtained within dissipative diabatic 
dynamics show the applicability of the theory to relatively light systems (Са + Ca) 
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as well as to medium (*Zr + ??Zr) [7| and heavy systems (РЬ + '?*Sn) [8]. In all cases 
the elastoplastic features of nuclei are easy seen. Within the framework of adiabatic 
single-particle motion the Far of motion for a single collective variable q reads: 


18B .. ду 
(00-1709 + dr: "d b (1) 


where the mass parameter B(q), the friction coefficient £(q) and the adiabatic poten- 
tial are defined within the shell model. However, the adiabatic approximation for the 
single particle motion is restricted to very small collective velocities [9,10]. A diabatic 
approximation which is defined by scaling the wave functions according to an irrota- 
tional flow imposed by the time dependence of the nuclear surface [1,5,11] has been 
find more realistic already for kinetic energies larger than 0.05-0.1 MeV per nucleon. 
This diabatic approximation to the single-particle motion in the time-dependent mean 
field is supported by time-dependent Hartree-Fock calculations [12]. On the basis of the 
general diabatic single-particle motion a transport theory has been formulated which is 
referred to as dissipative diabatic dynamics [3]. DDD results in the following changes 
of the collective equation of motion as compared to the adiabatic case (eq.(1)) : 

1. The strongly fluctuating adiabatic mass parameter is replaced by its smooth irro- 

tational value (B >В rrot). 


2. The adiabatic potential with shell corrections is replaced by its value at tempera- 
ture which results from excitation energy. This potential for temperatures largef 
than 2-3 MeV becomes close to the liquid drop potential. 


3. The Markovian friction is replaced by a retarded friction force 
e«t) — Јака). j 


The non-locality of the internal kernel K results from the equilibration within the in- 
trinsic degrees of fredom. By using a simple relaxation ansatz K becomes 


y= (t-t) 
K(t,t) = С exp EM, (3) 
where Tine (t) is the equilibration 370 This time has been estimated [13] within the 
Fermi-gas model to be 


Tinte(t) = t - 10-233 MeV/e'(t) , (4) 


where e'(t) is the time dependent excitation energy per nucleon and t*=20. The coef- 
ficient C(q) (a tensor in the case of more than one collective degree of freedom) stands 
for the stiffness of the system. 


The non-locality of the retarded friction term leads to the elastoplastic proper- 
ties of nuclei. If w denotes the frequency of the isoscalar vibrations (w = yc/ B) then 
for Tintr > WT! the system shows an elastic behaviour and for Tin & w”! the frictional 
limit is achieved with Markovian friction force (€ = C tintr)- 


The dissipative diabatic dynamics consists of two basic elements. Diabatic 
single-particle motion approximately describes the coherent quantum-mechanical cou- 
pling between intrinsic and collective degrees of freedom. The diabatic excitations 
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produce a repulsive force on the collective motion (the collective kinetic energy is pri- 
marily stored in the conservative potential). Then, dissipative collisions, essentially due 
to the residual two-body interactions, try to establish a new equilibrium distribution 
and destroy the diabatic potential. The intrinsic equilibration by the collisions is a 
time-irreversible process which leads to dissipation. It has been shown in [11] that this 
approach is applicable to nucleus-nucleus collisions in the range of 0.3 to 3 MeV per 
nucleon above the Culomb barrier. 


The aim of this paper is to discuss the problem of quasi-elastic recoil during 
central nucleus-nucleus collisions within one-dimensional as well as two-dimensional 
dynamics. The system 208РЬ + '*4Sn is chosen to maximize the effect of the elastic re- 
sponse [8]. In this case the adiabatic potential is more or less repulsive for all reasonable 
values of the collective variables and, hence, the elastic response cannot be obscured by 
the equilibrium forces. In cases of lighter systems ("Ca + “Са and Zr + Zr |7]) 
adiabatic potentials differ significantly. The minimum for the spherical configuration 
of the total system exists which produces a driving force towards fusion. For heavy 
systems the repulsive Coulomb force predominates and fusion is impossible. 


The paper is organized as follows. In sect. 2 the basic definitions and equations 
for a set of q = {qn} of collective variables are summarized. In sect. 3 the evaluation 
of the ingredients of the collective equation of motion is described. The results of one- 
dimensional trajectory calculations are discussed in sect. 4. In sect. 5 the preliminary 
results for two-dimensional dynamics of the same system are presented. 


2 Basic relations 


For a set q = {gn} of collective variables diabatic single-particle states are 
defined [5] by 


pa = exp |- 5C dent) Ми, | olesia) (5) 


where €, denotes the single-particle energies. All stationary states фо scale according 
to the same collective velocity field v = VW like: 


Fag | pa > = dl) V eV (Vw Ipa > (6 


for W = У, dnwa. The velocity field is allowed to describe compressions (AW # 0). 
With the scaling condition (6) all couplings proportional to g, vanish in the single- 
particle Schrodinger equation for the diabatic representation. Diabatic wave functions 
with different nodal structure (different character) exhibit crossings of the corresponding 
diabatic energy levels. The word 'diabatic! is derived from the greek word ‘тобабі сец’ 
for ‘to cross’. 


In the derivation of the transport equations [3] a time-smoothing procedure 
has been applied for the collective quantities which by definition should be the slow 
modes of the system. Within this formalism no ordinary friction term arises. Instaed, 
dissipation is obtained only through the changes of the occupation probabilities n, 
for the diabatic single-particle states. These changes are essentially due to two-body 
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collisions although one-body collisions from the remaining one-body coupling within 
the diabatic representation may also contribute. The collision term conserves the total 
mass and charge as well as the total energy. In practice we approximate the collision 
term by the relaxation equation 


dno(t) — (na(t) — %o(a, Т) m 
di A finie (E) ; : в 


where the intrinsic equilibration time (4) is time dependent via the changing excitation 
energy. The equilibrium occupation numbers Ж, are Fermi functions 


fa = {1 + exp [ee] | . | (8) 


The chemical potential u and the temperature T' are determined from the conservation 
of particle number and energy, respectively. 


The expectation value of the many-body hamiltonian, including two-body in- 
teractions, is given as sum of potential and kinetic collective energies 


1 oad 
< H >= V(q) a 2 у Anm . (9) 


The collective equations of motion (Euler-Lagrange equations) resulting from (9) are 


LOB aw. . 
aL Bać») — 2 € = Ё. = 72 Rol) ga. П (10) 


where the derivatives of Ње diabatic single-particie епегріез enter. Equations (7) and 
(10) form a set of coupled equations which are łocal in time, and hence Markovian. 
However, if we elimiaate the intrinsic variables п.(#) by formal integration of eq. (7) 


«(t = A) — [a exp |- J" oz] Sante) [ous (11) 
we obtain the non-Markovian collective equations of motion 


SE Bani) - z zz dnie Е 


for the collective variables. In the harmonic approximation the integral kernel is given 
by 
t 
Kot) = Com(a) exp |- 48 rzi,(6)] (13) 
and is referred to as the elastoplasticity tensor. The stiffness tensor is defined by 
С ae JA 
Caen) 09 


дд (ea 


and thus does not depend explicitely on t. We can easily recognize the elastic limit 
for пы, № (t — Ё) when K, (f, Г) becomes the stiffness tensor Cam = K,m(t,t' = t). 
For small amplitudes the corresponding vibrations can be identified wìth isoscalar giant 
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vibrations (1,3,4,14,15]. For słow motion where g(t) © q,,(t) we find the frictional 
limit with the friction tensor Enm = f dt Kam(t, Ё). 


The equilibrium force 


de 
K= Уһ) 5 . Qs) 
may be approximated by the derivative of the adiabatic potential which is smoothed 
like with a finite temperature. The corresponding potential should be very close to the 
liquid drop energy (the sum of surface plus Coulomb energies). 


3 Equilibrium force, mass and stiffness tensors for 
208pp + "Sp, 


In order to study the dynamic evolution of the colliding system according to (12) 
we have to calculate the ingredients of this equations. The basic problem of choosing 
the proper collective variables, however, arises. The most important collective variable 
is a quantity describing the relative motion of two colliding nuclei. As we deal micro- 
scopically with two-center shell model it is convenient to choose this quantity as ¢ — the 
distance between the centers of the two potentials (one can also relate it to the distance 
between the centres of masses). In the case of one-dimensional dynamics we merely 
consider the evolution of this single collective variable. 


Мы (© (MeV) 
Вх Ибо?) 


Fig.l. The adiabatic potential (solid line) and irrotational mass parameter 
(dot-dashed line) as functions of the elongation coordinate for ?*Pb + '?*Sn. 


Within the two-center shell model (TCSM) there also exist the other candidates 
for the collective variables. These are: deformation 6 = w,/w, (the same for both 
fragments), mass asymmetry and neck parameters — 6,a and e, respectively (for the 
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exact definitions see [5] ). Among them the deformation seems to be the most important 
and it is chosen to be a dynamical variable in the case of two-dimensional dynamics. 
The other parameters are fixed in calculations described below. 


The equilibrium force (15) is approximated by the derivative of the adiabatic 
potential with respect to the collective variables 


LL aV, 
PN ; hu--F н (16) 


The adiabatic potentiał is całculated as the sum of surface and Coulomb energies within 
a generalized liquid drop model with the shape defined by ¢ and 6. The calculations 
are performed according to the method of Moller and Nix [17]. In Fig.l the adiabatic 
potential for the system ?*Pb--!?*Sn is illustrated for the one-dimensional case (solid 
line). It is seen that apart from the narrow region of ¢ (slightly smaller than the distance 
at contact) the potential produces a repulsive force. The two-dimensional contour plot 
of the adiabatic potential for the same system is shown in Fig.2. 


Vadap(£,6) (MeV) 


0 5 10 15 
C (fm) 


Fig.2. Two-dimensional contour plot of the adiabatic potential for the sys- 
tem SPb + '?2*Sn. The energy separation between the dotted lines is 5 MeV. The 
thick dashed line shows the contact configurations for different shapes (sharp surface 
assumed). 


Note that 6 = 1,¢ = 0 corresponds to the spherical configuration of the total 
system, whereas 6 = 1 and large | ¢ | corresponds to separated spherical fragments. 
The dashed lines show the contact configuration for different deformations. The one- 
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dimensional case shown in Fig.1 corresponds to the cut of the Vi4(¢, ô) for 6 = 1. 


Bg 
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Br 
А2 И fr) 


Fig.3. Contour plots of the components of the mass tensor for the system 
208РЬ + '4Sn. In the upper part cuts are separated by 1 unit and the thick dashed 
line coincides with the cut B,, = 11. For the off-diagonal component the separation 
between dotted lines is 10 units and the one between solid lines is 40 units. For By, the 
separations are : 100 for dashed lines and 1000 for solid lines. The thick dashed line 


indicates the contact configurations for different shapes. 
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The mass tensor has already been studied in ref. [5] for various systems and it 
has been shown there that the mass tensor calculated from the cranking formula in the 
diabatic two-centre shell model is within 10% equal to the value given by the Werner- 
Wheeler approximation [16] to the irrotational flow (the cranking mass tensor which 
results from diabatic states exactly fulfilling the scaling condition (6) should coincide 
with its irrotational value). Therefore, in order to simplify time consuming calculations 
the Werner- Wheeler approximation has been used. In Fig.1 the diagonal B,, component 
of the mass tensor is also presented as a function of the collective variable ¢ (dot-dashed 
line). The mass parameter B,, is largest for separate nuclei (having the reduced mass 
value) and smoothly goes down for more compact configurations. For the same distance 
it is also smaller for prolate deformations (6j1) than for oblate ones (571). 


Fig.3 presents contour plots of all components of the mass tensor for the system 
under consideration. The dot-dashed line in Fig.l corresponds to the cut of the Bee 
according to 6=1. 
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Fig.4. The diagonal component of the stiffness tensor calculated for different 
values of the nuclear temperature. The solid line represents the results for T=2...3 MeV 
(deviations within linewidth). 


The stiffness parameter (calculated within DTCSM according to eq. (14)) as 
function of the distance between the centers of the single-particle potential is shown 
in Fig.4 (one-dimensional case) for different values of the nuclear temperature. The 
temperature is calculated via 


Т = \/8e*/MeV MeV (17) 


from the total excitation energy e* given by 


oo Bou - Vat) — Eus (t) 
A 
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where Ёсм, Ekin and A denote the bombarding energy in the centre-of-mass system, 
the collective kinetic energy and the total nucleon number, respectively. The stiffness 
parameter presented in Fig.4 exhibits strong shell effects for low temperatures which 
are smoothed out with increasing temperatures and practically vanish for T>2 MeV. 


Cg 
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Fig.5. Contour plots of the components of the stiffness tensor for the system 
208РЪ + !24Sn. In the upper part cuts are separated by 4 units between solid as well as 
dotted lines. For the off-diagonal component the separation between solid lines is 100 


units. In the lowest part the separations are : 1000 for dashed lines and 5000 for solid 
lines. The thick dotted lines indicate contact shapes. 
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In the two-dimensional case the components of the stiffness tensor are presented 
in Fig.5 for the case T=2 MeV when most of the shell effects vanish because for the 
energies under consideration the excitation energy is large enough (except perhaps from 
the very early stages of the collision process). One should notice that the non-diagonal 
elements of the stiffness tensor as well as of the mass tensor are not small and produce 
a strong coupling of the motion in the two collective variables. 


4 One-dimensional dynamics 


We have studied the dynamical evolution of the system 208 Pb+!**Sn as function 
of the incident energy above the Coulomb barrier, 


‚ _ Ecm — Veout 
in 7 Ara (19) 


between 0.1 and 5 MeV/u. 
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Fig.6. Quantities describing the time evolution of the system for three values 
of the incident energy above the Coulomb barrier, &;„ = 1 (dotted lines), 2 (dot-dashed 
lines) and 3 MeV/u (solid lines). Top: elongation ¢ as functin of time. Middle: intrinsic 
equilibration time Tin. Bottom: kinetic energy with respect to initial value. 
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Fig.6 shows the dynamicał evolution of the system for three values of the inci- 
dent kinetic energy 1, 2 and 3 MeV/u. For e,, =} and 2 MeV/u we get a remarkable 
elastic or quasi-ełastic behaviour of the system. In these cases the kinetic energies at the 
exit point (the barrier) are 50 MeV and 30 MeV, respectively, corresponding to ~65% 
and 225% of the initial energies. The interaction times (i.e. the times the system stays 
inside the Coulomb barrier) are almost the same (= 90. 10-235). The collisions proceed 
almost in the same way apart from the fact that for higher incident energy the inter- 
penetration of the colliding nuclei is deeper, while the excitation energy is higher. Since 
Tine i$ inversely proportional to e” it becomes smaller at higher bombarding energies. 
The qualitative behaviour of the trajectories change drastically as one goes to Ein = 3 
MeV/u. In this case the interaction time increases to more then three times the values 
obtained for e,, = 1 and 2 MeV/u. The reason is that the system is reflected back once 
by the elastoplastic force which is built up on the way towards reseparation. For times 
t>100:10-*5s practically all kinetic energy is dissipated and the nuclear temperature 
T>2 MeV is stabilized while the system creeps towards the Coulomb barrier. 


interaction tme 


t'an MeV/u) 


t'a (MeV/u) 


Fig.7. Kinetic energy at the exit point (left) and interaction time (right) as 
functions of the initial kinetic energy per nucleon above the Coulomb barrier for three 


different values of t*. The straight dotted line represents the initial kinetic energy (in 
the CM system). 


The results obtained for different values of the parameter t" are summarized in 
Fig.7 where the kinetic energies at the exit point (Coulomb barrier) and the interaction 
times are shown as functions of the incident kinetic energy for three different values of 
t*. For the cases t*=20 and é*=10 (the latter is close to the frictional limit) we have 
found a sudden change of the interaction time as function of the initial kinetic energy. 
There exists a particular initial energy (the value depending on t*) beyond which the 
interaction time becomes roughly three times larger and the quasi-elastic behaviour 
of the system switches to a frictional one where the kinetic energy at the exit point 
becomes very small. 
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5 Two-dimensional dynamics 


We consider here the central collisions of 7*Pb + '?*Sn at several incident 
energies in the range of 0.2 - 5 MeV above the Coulomb barrier. The time evolution 
of the system is obtained by solving the equation of motion (12) for two collective 
coordinates: elongation ¢ and deformation 6. The evolution of the system is calculated 
starting from the contact configuration of two spheres (6=1, ‹=12.6 fm) and assuming 
an initial velocity equal to the velocity of the relative motion (6(t=0)=0). Numerical 
integration of the equations of motion is performed using the predictor-corrector method 
and checking its accuracy by total energy conservation condition. 


5 10 ' 
$ (fm) C (fm) 


Fig.8. Two-dimensional trajectories for "Pb + !?*Sn for incident energies 
between 0.2 and 4 MeV/u above the Coulomb barrier. The solid line denotes the 
trajectory inside the barrier (compound configuration) and the dotted line shows the 
evolution after reseparation. The dashed line indicates contact configurations. 


The trajectories of the colliding system for different incident energies are illus- 
trated in Fig.8. All results correspond to t*=20. The interpenetration of nuclei increases 
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with the increasing energy and for £i > 4 MeV/u ç even becomes negative for some 
time. 


The shapes of the trajectories are easy to understand if one looks at the contour 
plots of V.u(¢,6) and С;‹(,5). It is seen that the system evolves along the valley of the 
smallest stiffness along which the dynamically built up elastoplastic forces are smallest. 
Along this valley the equilibrium forces are also small (see Fig.2). The coordinates, 
however, are strongly coupled by the non-diagonal elements of the tiffness and mass 
tensors. For low energies this heavy system exhibits similar quadrupole vibrations as 
the lighter ones (Ca + “Ca and ®Zr + Zr [7]). For energies ei, > 4 MeV/u, our 
investigation is not yet finished because as can be seen from the lowest-right part of Fig.8 
the system escapes from the region where the ingredients of the equations of motion 
have been calculated microscopically (however is still inside the barrier). 


Energy 


— t'=10 
— t20 
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Fig.9. The sum of the kinetic and potential energies (at the exit point at the 
barrier) with respect to its initial value as function of the incident energy for different 
values of the parameter t”. 


As the Fermi-gas model estimate of the intrinsic equilibration time is not well 
determined we performed additional calculations with t*—10. The final results are 
collected in Figs.9 and 10. These results qualitatively agree with those of Fig.7 obtained 
with one-dimensional dynamics. Fig.9 shows that for lower incident energies less than 
one half of the total energy is dissipated whereas for higher incident energies almost all 
energy is dissipated. 


The interaction time shown in Fig.10 is almost the same for lower energies and 
increases rapidly for higher ones. This result is also in qualitative agreement with the 
one of one-dimensional dynamics. Quantitative results and the whole picture of the 
collision process within two-dimensional dynamics, however, differ much. 
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Fig.10. The interaction time (the time the system stays iaside the Coulomb 
barrier) as function of the incident energy. 


In closing we note that our studies still cannot give quantitative predictions 
because we need a wider space to calculate the trajectories for higher incident ener- 
gies. The other collective degrees of freedom (mass asymmetry, neck) might also be 
important. 
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Streszczenie 


W pracy prezentowane są obliczenia nuneryczne dotyczące 
centralnych zderzeń jąder z dyssypacją energii. Głównym celem 
pracy jest przedyskutowanie problemu kwazielastycznego odrzutu 
z wykorzystaniem jednowymiarowej i dwuwymiarowej dynamiki, 
Ewolucję zderzających się jąder opisano przez rozwiązanie kla- 
sycznych równań ruchu dla współrzędnych kolektywnych, Jako 
współrzędne kolektywne przyjęto zmienną opisującą ruch względny 
fragmentów oraz deformację kwadrupolową, Tensor masowy 1 siły 
równowagi obliczono makroskopowo /uogólniony model kropli cie- 
czy/, podczas gdy elastoplastyczność wzięto pod uwagę w modelu 
mikroskopowym /diabatyczny, dwucentrowy model powrokowy/. 
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Rezultaty dynamiki jednowymiaroxej jakościowo zgodne są 
z rezultatami otrzymanymi w przypadku dwuwymiarowym. Główną 
konkluzja pracy jest wniosek, że w centralnych zderzeniach 
ciężkich jąder dla niskich energii zderzeń /1-2 MeV/u/ można 
oczekiwać efektów kvazielastycznych. 


РЕЗЮНЕ 


В рабсте призодятся численные расчеты связанные с централь- 
ными столкновениями ядер с диссипацией энергии. Глазной целью 
работы является обсухдение проблемы квазиэластичной отдачи с HC- 
пользованием одно- и двухмерной динамик. Эволюция столкивающихся 
ядер описывается путем решения классических уравнений движения 
для коллективных координат. В качестве коллективных координат 
приняты переменная списызающая относительное двизение фрагментов 
и квадрупольная дефорыеция. массовый тензор и силы равновесия 
вычислены хакроскопически (обобщенная модель жидкой капли), в TO 
время как эластопластичность принималась во внимание в микроско- 
пической модели (адиабатическая двух-центричная оболочечная мо- 
дель). 

Результаты одномерной динамики согласуются качественно с 
данными полученными в двухмерном случае. Основным заключением 
работы является то, что в центральных столкновениях тяжелых 
ядер при низких энергиях столкновений (I-2 МэВ/н) можно ожидать 
квазиэластичных эффектов. 


